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Abstract. Absolute photoionization cross-section calculations are presented for 
Se"*" using large-scale close-coupling calculations within the Breit-Pauli and Dirac- 
Coulomb R-matrix approximations. The results from our theoretical work are 
compared with recent measurements PQ .2, \3\ made at the Advanced Light Source 
(ALS) radiation facility in Berkeley, California, USA. We report on results for 
the photon energy range 18.0 eV - 31.0 eV, which spans the ionization thresholds 
of the 4 S2,„ ground state and the low-lying 2 D°,„ ,,„ and 2 P° ,„ . ,„ metastable 
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states. Metastable fractions are inferred from our present work. Resonance ener- 
gies and quantum defects of the prominent Rydberg resonances series identified 
$-H in the spectra are compared for the 4p — > nd transitions with the recent ALS 

Vyi experimental measurements made on this complex trans-iron element. 
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1. Introduction 

Recently, detailed photoionization cross-section measurements for the trans-iron 
element Se + , carried out at the Advanced Light Source (ALS) synchrotron radiation 
facility in Berkeley California were reported in the literature [TJ [5J |3| . The photon 
energy studied in the experiments ranged from 18.0 eV to 31.0 eV and was scanned with 
a nominal spectral resolution of 28 meV. Measurements were also made at a higher 
photon energy resolution of 5.5 meV from 17.75 to 21.85 eV spanning the 4s 2 4p 3 4 Sg , 2 
ground-state ionization threshold and the 2 P§/2, 2 P1i 2 , 2 Dg/ 2 , and 2 D°, 2 metastable 
state thresholds. Although detailed resonance structure in the photoionization spectra 
was reported, it was not possible to positively determine the metastable components 
in the ion beam. Photoionization cross section calculations have been made recently 
on this complex ion using AUTOSTRUCTURE g] within the Multi-Configuration- 
Breit-Pauli (MCBP) distorted-wave approximation for a wider range of energies and 
for higher charged states. It was concluded from that study on Se ions [I] (by visual 
comparison) that best agreement with experiment was obtained by weighting the 
contribution of the ground configuration states ( 4 Sg\ 2 , 2 Dg/ 2 , 2 Dg, 2 , 2 P1/ 2 , 2 P-j/ 2 ) by 
(0.53, 0.15, 0.05, 0.11, 0.16) respectively. As pointed out by Sterling and Witthoeft 
U, R matrix techniques are necessary to reproduce resonance interference effects. In 
the present work we report on large scale photoionization cross section calculations 
for this complex trans-iron ion in both the Breit-Pauli and Dirac Coulomb R-matrix 
approximations. We infer metastable fractions and benchmark our theoretical results 
with the available experimental data [TJ El |3] providing confidence in the data on this 
trans-iron element for applications [SI [5] . 

Photoionization (PI) of atomic ions is an important process in determining the 
ionization balance and hence the abundances of elements in photoionized astrophysical 
nebulae. It has recently become possible to detect neutron(n)-capture elements 
(atomic number Z > 30, e.g. Se, Kr, Br, Xe, Rb, Ba and Pb) in a large number of 
ionized nebulae j7j [8j [9]. These elements are produced by slow or rapid n-capture 
nucleosynthesis (the "s-process" and "r-process," respectively) [SJ QJj]. Measuring 
the abundances of these elements helps to reveal their dominant production sites 
in the Universe, as well as details of stellar structure, mixing and nucleosynthesis 
1 11 L H2l fT3l [T4l [TBI [TBI IT71 [5]. These astrophysical observations are the motivation to 
determine the photoionization and recombination properties of n-capture elements. 

Various n-capture elements have been detected in the spectra of planetary nebulae 
(PN^p, the photoionized ejecta of evolved low- and intermediate-mass stars (1-8 solar 
masses) . Planetary Nebula progenitor stars may experience s-process nucleosynthesis 
|141 118|. in which case their nebulae will exhibit enhanced abundances of trans-iron 
elements. The level of s-process enrichment for individual elements is strongly sensitive 
to the physical conditions in the stellar interior [T4l [T6] . 

The underlying difficulty in studying s-process enrichment in PNe is due to the 
large uncertainties (factors of 2-3) of n-capture element abundances derived from 
the observational data [2J |SJ. There are two main causes for these uncertainties. 
First, the low cosmic abundances of trans-iron elements, means their emission lines 
are sufficiently weak that generally only one or two ions of each element can be 
detected in individual PNe. Second, while robust ionization corrections can be derived 

§ Following astrophysical nomenclature, we abbreviate the singular form planetary nebula as "PN" 
and the plural form as "PNe." 
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from numerical simulations of nebulae [19], [20]. this method relies on the availability 
of accurate atomic data for processes that affect the ionization equilibrium of each 
clement. In photoionized nebulae, these atomic data include photoionization cross- 
sections and rate coefficients for radiative and dielectronic recombination and charge 
exchange reactions are unknown for nearly all n-capture element ions. Uncertainties 
in the photoionization and recombination data of n-capture element ions can result in 
elemental abundance uncertainties of a factor of two or more [7]- 

The present theoretical investigation is a study to determine the photoionization 
and recombination properties of n-capture element ions, motivated by the 
astrophysical detection of these species and the importance of their elemental 
abundances for testing theories of nucleosynthesis and stellar structure. Determining 
the appropriate atomic data over the range of energies and temperatures 
encountered in astrophysical environs necessitates a predominantly theoretical 
approach. Experimental measurements are helpful to validate theoretical studies, 
particularly in the case of complex heavy ionic systems such as low charge states 
of trans-iron elements. We note that experimental studies are often limited in the 
photon energy range studied, energy resolution and to the charge states able to be 
investigated. In addition, many of the experiments carried out are contaminated by 
metastable states and require theoretical estimates to help unravel the components in 
the ion beam. 

The prime motivation for choosing Se as the first element for investigation is 
due to the fact it has been detected in nearly twice as many PNe (70 in total) as 
any other trans-iron elements (9] in addition to recent experimental measurements 
becoming available to compare with [TJ |5] . Experimental photoionization studies of 
other astrophysically observed n-capturc elements have already been conducted by 
various groups in the cases of select Kr [211 |2"21 123"] and Xe ions [241 125] , 

In this paper we present theoretical determinations of the absolute Se + 
photoionization cross-section near the ground-state ionization threshold and compare 
our results with recent experimental measurements [TJ [5] for this trans-iron ion. We 
note that the results from our work can be applied to calibrate a broader theoretical 
effort to determine the photoionization and recombination properties of n-capture 
element ions over a large range of energies and temperatures. Our data would be 
suitable for incorporation within photoionization modelling codes such as Cloudy and 
XSTAR[T!?l I2TJ] . which numerically simulate the thermal and ionization structure of 
astrophysical nebulae. Photoionization models can be used to derive accurate and 
robust ionization corrections for trans-iron elements, thereby enabling the accurate 
determination of their abundances in the Universe. We also intend to explore the 
properties of other members of this iso-electronic sequence, for example, Br 2+ and 
Rb 3+ along with other low-charged Se ions in subsequent publications. 

The layout of this paper is as follows. Section 2 presents a brief outline of 
the theoretical work. Section 3 presents our theoretical PI cross sections results 
where comparisons are made with the recent Se + photoionization experimental 
measurements [U [2 |3]- Resonance energies and quantum defects for the 
various dominant Rydberg series (seen in both the experimental and theoretical 
photoionization spectra) are tabulated and compared. Finally in section 4 conclusions 
are drawn from the present investigation. 
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2. Theory 

Relativistic corrections can play an increasingly significant role for atoms with charge 
Z > 18 [25] . It was shown in recent work on electron-impact excitation of neutral neon 
(Z = 10) that a jK coupling scheme was required to accurately describe even the low 
lying energy target levels [53]. The semi-relativistic Breit-Pauli R-matrix approach 
has been successfully applied repeatably to atomic ions below Z=30 for numerous 
atomic systems by the Iron Project [27]. For atomic ions with Z > 30, a Dirac- 
Coulomb approach to electron and photon interactions is ideal, due to the importance 
of relativistic and correlation effects in accurately modelling these systems 

Computationally, the increased number of relativistic orbitals over the 
corresponding non-relativistic orbitals can cause a dramatic increase both in the 
number of integrals required and in the complexity of the angular algebra. We note 
that in a recent work for the electron-impact excitation of Fe 4+ [25] and Fe 5+ [2S] 
hundreds of millions of Racah coefficients were calculated for the relativistic treatment 
of this open-d shell system. To address the future challenge of electron scattering in 
which hundreds of levels and thousands of scattering channels are to be considered an 
efficient parallel version [30J of the DARC [3"T ] l3"2" | l3"3" ] 13"4"] suite of codes was developed. 

A natural extension of this work was the formation of the numerous bound-free 
dipole matrix elements required for photoionization in an equally efficient manner, 
distributed over an equivalent number of processors. This enables photoionization 
and photo-recombination calculations on complex halogen- like ions for example, Ar + , 
Kr + and Xe + , to be carried out at the same degree of accuracy [33] as those for 
electron impact excitation. It should be noted that in addition to large scale relativistic 
photoionization calculations, the same dipole matrix elements are the foundation of 
radiationally damped electron-impact excitation studies required in highly charged 
ions. With these extensions to the DARC codes one can address the complex problem 
of trans-iron element single photon ionization. The flowchart in figure [II with brief 
annotations beside each module of the suite of DARC codes, illustrates the R-matrix 
sequence required to obtain photoionization cross sections with the parallel Dirac 
Atomic R-matrix Codes. These recent modifications now include the concurrent 
formation of every dipole matrix, reducing the total time substantially required for 
the formation of all dipoles to that required by a single dipole. These developments 
open the frontier to enabling comprehensive pioneering large scale calculations along 
iso- nuclear sequences to be performed now within a feasible time frame. 

Using this suite of DARC codes detailed photoionization cross section calculations 
for the Se + ion were performed for the ground and all the excited metastable levels 
associated with the 4s 2 4p 3 configuration. To benchmark our theoretical results we 
compare them with recent high resolution experimental measurements made at the 
Advanced Light Source synchrotron radiation facility [U E] [3] . 

In R-matrix theory, all photoionization cross section calculations require the 
generation of atomic orbitals based primarily on the atomic structure of the residual 
ion. The present theoretical work for the photoionization of the Se + ion employs 
relativistic atomic orbitals generated for the residual Se 2+ ion, which were calculated 
using the extended-optimal-level (EOL) procedure within the GRASP structure code 
[36} 137} 125] . We explored various structure and scattering models of increasing 
complexity with the R-matrix approach in order to investigate an optimum model 
for determining the total photoionization cross sections. 

Photoionization cross section calculations for this complex trans-iron system 
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Figure 1. Flowchart illustrating the steps necessary to obtain photoionization 
cross sections with the parallel Dirac Atomic R-matrix Codes (DARC). The codes 
are available at http://connorb.freeshell.org/codes/DARCdipole/ 



began by including only levels from the 4s 2 4p 2 , 4s4p 3 and 4p 4 configurations which 
gave rise to 20 Jit levels. We found this model to be inadequate to represent the 
complex experimental spectrum and it was necessary to extend the theoretical model. 
Our first approach was to supplement the configurations of this 20 level model with 
addition levels from several 3d-hole configurations, 3d 9 4s 2 4p 3 , 3d 9 4s4p 4 and 3d 9 4p 5 . 
This increased the number of levels in our close coupling expansion to 126 yielding 
improved agreement with experiment. Further extensive theoretical models were then 
investigated with the inclusion of levels from configurations involving a 4d orbital 
in our close coupling expansion rather than levels from the 3d hole configurations. 
The first of these models included all 246 levels arising from the eight configurations: 
4s 2 4p 2 , 4s4p 3 , 4p 4 , 4s 2 4p4d, 4s 2 4d 2 , 4p 2 4d 2 , 4s4p 2 4d and 4p 3 4d were included in 
the close-coupling expansion. The second model investigated was when this 246-level 
model was augmented with levels from an additional 4s4p4d 2 configuration, giving 
rise to a total of 336 levels. 

The final model augmented this 336-level model with the 4s 2 4p5s, 4s 2 4p5p and 
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Figure 2. (Colour online) Photoionization of the 4 SS/ 2 ground state of the Se+ 
ion. The cross section from threshold to a photon energy of 30 eV are illustrated. 
Theoretical results have been obtained from both the DARC and Breit-Pauli 
approximations with an increasing number of levels included in the close-coupling 
approximations. The theoretical cross sections have been convoluted with a 28 
meV FWHM gaussian distribution for comparison purposes. 



4s 2 4p5d configurations which gave rise to a total of 362-levels. Table 1 shows a sample 
comparison of our results for the lowest 12 levels of the residual Se 2 (Se III) ion with 
the available experimental and theoretical work. 

Photoionization cross section calculations were performed in the Dirac-Coulomb 
and Breit-Pauli approximation with the various models. The advantage of the DARC 
calculation over the Breit-Pauli result is marginal for the photoionization of outer shell 
electrons. However, it does provide the repeatability of the result using completely 
independent structure codes and for the most part independent collisional codes to 
provide the cross section. 

Considering the ionization potential between the ground state of Se II and that of 
Se III is approximately 1.6 Ryds and according to NIST tables the 4s 2 4p4d, 4s 2 4p5^ 
levels are over 1 Ryd above the Se III ground state. This gives a photon energy 
of greater than 36 eV to photoexcite to either the 4d 2 4p4d or the 4s 2 4p5s levels. 
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Table 1. Comparison of the theoretical energies using GRASP and the Breit- 
Pauli approximation from the 362-level model with previous work, relative 
energies are in Rydbergs. A sample of the first 12 levels of the Se III ion are 
compared with AUTOSTRUCTURE g] and the NIST [38] tabulated values. 



Level 


Config. 


Term a 


Energy 


Energy 


Energy 


Energy 


A(%) b 








NIST 


GRASP 


Breit-Pauli d 


AUTOS 6 




1 


4s 2 4p 2 


3 Po 


0.000000 


0.000000 


0.00000 


0.00000 


0.0 


2 


4s 2 4p 2 


3 Pi 


0.015870 


0.015061 


0.01419 


0.01580 


5.1 


3 


4s 2 4p 2 


3 P 2 


0.035880 


0.035409 


0.03347 


0.03620 


1.3 


4 


4s 2 4p 2 


T 2 


0.118760 


0.140979 


0.13835 


0.14370 


18.7 


5 


4s 2 4p 2 


So 


0.259070 


0.285184 


0.27116 


0.25950 


10.1 


6 


4SMP 3 


5 S°, 


0.579910 c 


0.566575 


0.55304 


- 


2.3 


7 


4s 1 4p 3 


3 D? 


0.830080 


0.853381 


0.81498 


0.83080 


2.8 


8 


4s 1 4p 3 


3 D§ 


0.844950 


0.854258 


0.81543 


0.84500 


1.1 


9 


4SMP 3 


3 D§ 


0.879810 


0.859092 


0.81900 


0.87980 


2.3 


10 


4s 1 4p 3 


3tjo 
r 


0.970335 


0.991048 


0.94428 


0.96750 


2.1 


11 


4s 1 4p 3 


3T30 

r 2 


0.970636 


0.994054 


0.94614 


0.96540 


2.4 


12 


4s 1 4p 3 


3po 


0.971329 


0.992406 


0.94616 


0.96760 


2.2 



a 2S+1 tit 
J 

b GRASP, absolute percentage difference relative to NIST values 

c Esteves PhD thesis [J 

d Breit-Pauli 

Sterling and Witthoeft g] 



Therefore, their influence on the ground state photoionization cross section within the 
energy range presented in this paper, is solely in terms of changes in structure through 
CI mixing. 

Figure 2 shows that the inclusion of the additional 5s, 5p and 5d orbitals in 
the basis set has a minimal effect on the total photoionization cross section for the 
4 Sg , 2 ground state over the energy range investigated and for energies less than 31 
cV. However, these comparisons with experiment are slightly misleading as the energy 
levels of the DARC R-matrix calculation are shifted to experimental values prior to 
the matrix diagonalization. Since converged results are obtained within the 336 level 
model we used this approximation to determine all the remaining PI cross sections for 
the metastable states with the DARC codes. 

The R-matrix boundary radius of 8 Bohr radii was sufficient to envelop the radial 
extent of all the atomic orbitals of the residual Se 2+ ion. A basis of 16 continuum 
orbitals was sufficient to span the incident experimental photon energy range of up 
to 32 eV. The 336-state model produced a maximum of 1758 coupled channels in our 
scattering work with Hamiltonian matrices of dimension of the order of 28,300 by 
28,300 in size. Due to dipole selection rules, for total ground state photoionization we 
need only consider the three bound-free dipole matrices, 2J 7r = 3° — > 2J 71 " = l c ,3 e ,5 e 
whereas for the excited metastable states one requires 2J W = 5° — >• 2J 7r = 3 C , 5°, 7 e the 
2J 7r = 3° ->• 2J ff = l c , 3°, 5° and 2J 7r = 1° -> 2J W = 1°, 3° bound-free dipole matrices. 
For each initial state the outer region electron-ion collision problem was solved (in 
the resonance region below and between all thresholds) using a suitably chosen fine 
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energy mesh of 5xl0 -7 Rydbergs (« 6.8 ^eV) to fully resolve all the extremely fine 
resonance structure in the appropriate photoionization cross sections. The multi- 
channel R-matrix eigenphase derivative (QB) technique (applicable to atomic and 
molecular complexes) of Berrington and co-workers [39, 40 41J was used to determine 
the resonance parameters. In the QB method, the resonance width V is determined 
from the inverse of the energy derivative of the eigenphase sum S at the position of 
the resonance energy e r via 

-l 



r = 2 



dS 



dE 



= 2[<J']iL . (i) 



E 
E=e r 



Prominent photoexcitation-autoionization Rydberg resonances series are analyzed and 
identified as 4p — > nd transitions, and the metastable fractions of the experimental 
Se + ion beam arc estimated. Finally, in order to compare directly with recent 
experimental work on this system [U E], the theoretical cross sections were 
convoluted with a gaussian function of appropriate width (FWHM of 28 meV and 5.5 
meV respectively) and weighted accordingly in order to simulate the energy resolution 
of the measurements performed at the Advanced Light Source (ALS) synchrotron 
radiation facility. 

3. Results and discussion 



Figure [2] shows the convergence of the PI results obtained from the Breit-Pauli and 

5 3/2 



Dirac-Coulomb (DARC) approximations for the 4s 2 4p 3 4 S2 , 2 ground state using an 



n=4 basis. The theoretical PI cross section results are illustrated for the 4s 2 4p 3 4 Sg , 2 
ground state of this complex trans-iron element with an increasing number of levels 
retained in the close-coupling expansion. As illustrated in figure [2] the 336-state 
model yields a suitable description of these half collision processes for the Se + ion, 
where converged results are obtained with a larger close-coupling approximation. At 
this level of approximation it ensures consistency in our theoretical work. Using this 
336-level model we obtained the PI cross sections for all the remaining metastable 
levels associated with the 4s 2 4p 3 configuration within the Dirac-Coulomb (DARC) 
approximation. The theoretical PI cross section results for the individual levels arising 
from the 4s 2 4p 3 configuration over the photon energy range 18 eV to 31 eV from the 
Dirac-Coulomb (DARC) approximation are shown in figure ^ where the results are 
convoluted with a gaussian distribution of 28 meV FWHM in order to compare with 
the high resolution ALS experimental results. The ALS measurements for Se + [2] 
were taken over the 18 eV to 31 eV photon energy range at an energy resolution of 
28 meV, these measurements are illustrated in the bottom panel of figure [3] The 
absolute photoionization cross sections were measured at discrete photon energies and 
are shown in the figure by solid circles with their associated uncertainties. 

To determine the fractions of metastable ions present in the experimental ion 
beam photoionization cross section calculations were performed using the Dirac- 
Coulomb R-matrix method for both the 4 Sg , 2 ground state and the excited metastable 

states ( 2 Dg /2 , 2 D£ , 2 , 2 P° /2 and 2 P3/ 2 ) associated with the 4s 2 4p 3 configuration. The 
individual fractions were estimated by first scaling the highest-lying metastable state 
to the absolute experimental cross section. Each subsequent state was then scaled 
such that the sum of the theoretical cross sections matched the experiment as closely 
as possible, with all four metastable states together with the ground state fractions 
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Figure 3. (Colour online) Theoretical photoionization cross sections (336-level 
close-coupling DARC calculations) for the 4 S?, 2 ground state, 2 D?, 2 , 2 D°,„, 

2 P°,„ and 2 P?/ 2 metastable states arising from the 4s 2 4p 3 configuration are 
illustrated (Top 5 panels). The theoretical results were convoluted with a gaussian 
of 28 meV FWHM to simulate the photon energy resolution of the experimental 
measurements (bottom panel) taken at the Advanced Light Source synchrotron 
radiation facility. 



3 , 2 , -^ 3 / 2 , ^5/ 2 ; * 1/2 allu x 3/2 states by (0.52, 0.13, 0.19, 0.03 and 0.13) 



ultimately constrained to sum to unity. Using this technique the various fractions 
were determined and it was found that a non-statistical distribution of the ground and 
metastable states gave best agreement with experiment by weighting the contribution 
of the 4 S° /2 , 2 ~D°m, 2 D? /0 , 2 P? rt and 2 P° 
respectively. 

As noted previously, photoionization cross section calculations have been made 
recently on this complex ion using AUTOSTRUCTURE 0] within the Multi- 
Connguration-Breit-Pauli (MCBP) distorted-wave approximation for a wider range 
of energies and for higher charged states of the ion. Using this distorted wave 
approximation the authors find that the direct PI cross sections for these Se ions differ 
by 30% to 50% when compared to the available ALS experimental measurements. 
It was concluded from that study on Se ions H] (by visual comparison) that best 
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Figure 4. (Colour online) Comparison of ALS photoionization cross-section 
measurements from 18 eV to 21 eV taken at 28 meV (shaded area) with the 
336-state close-coupling DARC calculations. The theoretical cross sections (solid 
black line) for each of the five initial states have been convoluted with a 28 meV 
FWHM gaussian and a weighting of the states (see text for details) to simulate 
the measurements. The absolute measurements (solid black circles) have been 
obtained with a larger energy step. The error bars give the total uncertainty of 
the experimental data. The various Rydberg series are indicated 



agreement with experiment was obtained by weighting the contribution of the ground 
configuration states ( 4 S° /2 , 2 D° /2 , 2 D° /2 , 2 P° /2 , 2 P° /2 ) by (0.53, 0.15, 0.05, 0.11, 0.16) 
respectively. Therefore a non statistical distribution of the PI cross sections in both 
theoretical approaches is seen to provide the best agreement with experiment. 

Furthermore the Fano profile in the PI cross section near 23.0 eV is not reproduced 
by the AUTOSTRUCTURE calculations [4, (R matrix techniques are necessary to 
reproduce resonance interference effects) , and the lack of this interference causes the 
discrepancy in the direct cross section in this energy region. The calculated and 
experimental Se + direct cross sections agree to within 30% to 50% over the range 
of experimental energies, except above 30 eV, where the experimental measurements 
are more uncertain since only a single absolute measurement was made above 27 eV 
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Figure 5. (Colour online) Comparison of ALS photoionization cross-section 
measurements from 18 eV to 21.8 eV taken at 5.5 meV (shaded area) with the 
336-state close-coupling DARC calculations. The theoretical cross sections (solid 
black line) for each of the five initial states have been convoluted with a 5.5 meV 
FWHM gaussian and a weighting of the states (see text for details) to simulate 
the measurements. The absolute measurements (solid black circles) have been 
obtained with a larger energy step. The error bars give the total uncertainty of 
the experimental data. The various Rydberg series are indicated. 



(compared to absolute measurements at lower energies, where independent estimates 
were performed 2 to 3 times for each energy). Therefore, they report uncertainties 
in the experiment of over 50 % in the high energy region. As R matrix techniques 
are necessary to reproduce resonance interference effects we have carried out these 
calculations in the present investigation. 

In figure |4j the rich resonance structure observed in the experimental spectra is 
generally well represented in the metastable region by theory. This is still apparent 
when the theoretical work is compared with the experimental measurements taken at 
an even higher resolution [3] . Figure p] shows the ALS experimental measurements 
taken at a higher resolution of 5.5 meV (from 18 eV to 21.8 eV) together with the 
theoretical results convoluted with a gaussian of the same FWHM to simulate the 
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Figure 6. (Colour online) Comparison of ALS photoionization cross-section 
measurements from 21 eV to 31 eV taken at 28 meV (shaded area) with the 
336-state close-coupling DARC calculations. The theoretical cross sections (solid 
black line) for each of the five initial states have been convoluted with a 28 meV 
FWHM gaussian and a weighting of the states (see text for details) to simulate 
the measurements. The absolute measurements (solid black circles) have been 
obtained with a larger energy step. The error bars give the total uncertainty of 
the experimental data. 



experimental resolution and non-statistically weighted as before. Here again there is 
good agreement between the experimental and theoretical results. 

In the energy region 21 eV to 31 eV, (figure [6]) the theoretical results are larger 
than experiment. The solid line corresponds to the sum of the ground and metastable 
state fractions of the theoretical calculations. As indicated in the ALS experimental 
studies for Se + [TJ [2], contamination from higher order radiation produced by the 
undulator beamline may be a possible cause for uncertainties. Previous experiments 
on Xe 3+ [32] estimated a lower fraction of 2 % at around 40 eV. The measurements on 
Se + [TJ [2] are between these two photon energies and it was concluded in that work 
[TJ [5] the maximum contamination is no greater than 4%. At lower photon energies, 
the contamination of higher order radiation is expected to be larger, but not by more 
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Figure 7. (Colour online) Eigenphase sum <5 and its derivative 5' (for each 
2J 7r symmetry contributing to the PI cross section a for the 4s 2 4p 3 2 P?/ 2 initial 

state) as a function of photon energy in the region below the *D 2 and So excited 
state thresholds of the residual Se 2 + ion for (a) 2J=1 C symmetry, (b) 2J=3 e 
symmetry and (c) PI cross section a for the 2 P?/ 2 metastable state convoluted 
with a gaussian of 28 meV FWHM (solid black line). The prominent Rydberg 
resonances series in the PI cross section have been identified (open and solid 
triangles) as members of two different Rydberg autoionizing series converging to 
the X D 2 and 1 Sq excited states of the residual Se 2 + ion indicated by thick vertical 
lines. 



than a factor of 2D 3 compared to the contamination at 40 eV. It was concluded from 
the Se + ALS measurements [TJ [5] that the experimental uncertainties on the absolute 
measurements made are estimated to be 30% including the possible contamination of 
the photon beam by higher order radiation. 

Earlier studies on the photoionization spectrum of + at the Advanced Light 
Source synchrotron radiation facility fl3~| were used to identify the resonance structure 
due to the metastable states jT] |2j. This is apparent when the spectra of + and 
Se + are compared, the resemblance is remarkable but not totally unexpected. These 
respective ions have an np 3 configuration (n=2 for + and n=4 for Se + ) in their 
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Figure 8. (Colour online) Eigenphase sum <5 and its derivative 8' (for each 2J' r 
symmetry contributing to the PI cross section a for the 4s 2 4p 3 2 P° ,„ initial state) 

as a function of photon energy in the region below the 1 D2 and 1 So excited state 
thresholds of the residual Se 2 + ion for (a) 2J=1 C symmetry, (b) 2J=3° symmetry, 
(c) 2J=5 C symmetry and (d) PI cross section a for the 2 P < i/ 2 metastable state 
convoluted with a gaussian of 28 meV FWHM (solid black line). The prominent 
Rydberg resonances series in the PI cross section have been identified (open and 
solid triangles) as members of two different Rydberg autoionizing series converging 
to the X D2 and 1 So excited states of residual Se 2 + ion indicated by thick vertical 
lines. 



outermost shell giving rise in each case to the same term designation for the 4 Sg , 2 



ground state and the 2 ^>%/ 2 5/2 an d 2 P°/2 3/2 



metastable states. In the case of O , the 
region of the spectrum covering the energy range from the 2 P°/ 2 3 / 2 threshold to the 
4 Sg j 2 threshold consisted of very strong resonances due to 2p — > nd electron excitations 
of the metastable ions as well as weak resonances due to 2p — > ns transitions. Similarly, 
the Se + spectra, shows the corresponding 4p — > nd transitions but the 4p — > ns series 
appear to be too weak to be observed. 

At energies below 21.2 eV there are prominent Rydberg resonance features 
resulting from photoionization of the initial metastable Se + states in the experimental 
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Figure 9. (Colour online) Eigenphase sum <5 and its derivative <$' (for each 2J' r 
symmetry contributing to the PI cross section a for the 4s 2 4p 3 2 D° ,„ initial state) 

as a function of photon energy in the region below the 1 D2 and 1 So excited state 
thresholds the residual Se 2 + ion for (a) 2J=1° symmetry, (b) 2J=3 e symmetry, 
(c) 2J=5 C symmetry and (d) PI cross section a for the 2 D° ,„ metastable state 
convoluted with a gaussian of 28 meV FWHM (solid black line). The prominent 
Rydberg resonances series in the PI cross section have been identified (solid 
triangles) as members of a Rydberg autoionizing series converging to the *D2 
threshold of the residual Se 2 + ion. The weak resonance series converging to the 
^o threshold is not observed in the PI cross sections. 



and theoretical cross sections, as is seen in figure [3] Multiple overlapping Rydberg 
resonance series in the spectra are also clearly visible converging to the 4s 2 4p 2 1 D2 and 
4s 2 4p 2 ^o thresholds respectively of the residual Se 2+ ion. The presence of interloping 
resonances disrupt the regular Rydberg pattern of resonances, as can be seen in the 
spectra for the 2 P°/ 2 and ^3/2 metastable states, respectively figures w\ and pi 

To obtain the resonance parameters a theoretical analysis with a niocuncation 
to the widely used eigenphase derivative (QB) technique (for jj -coupling) that 
determines resonance parameters in electron collisions with atomic and molecular 
systems [22] SB HI] was used. This technique exploits the properties of the R-matrix 
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Figure 10. (Colour online) Eigenphase sum 5 and its derivative 8' (for each 2J' r 
symmetry contributing to the PI cross section a for the 4s 2 4p 3 2 D° ,„ initial state) 

as a function of photon energy in the region below the 1 D2 and 1 So excited state 
thresholds of the residual Se 2 + ion for (a) 2 J=3 e symmetry, (b) 2J=5 e symmetry, 
(c) 2J=7 C symmetry and (d) PI cross section a for the 2 D° ,„ metastable state 
convoluted with a gaussian of 28 meV FWHM (solid black line). The prominent 
Rydberg resonances series in the PI cross section have been identified (solid 
triangles) as members of a Rydberg autoionizing series converging to the D2. 
The weak resonance series converging to the 1 So threshold is not observed in the 
PI cross sections. 



in multi-channel scattering and is an excellent method for locating and determining 
the properties of narrow resonances. For this trans-iron element our theoretical results 
(obtained for the dominant resonance features in the photoionization spectra) are 
compared with the recent ALS measurements [1, 2J. From the ALS measurements 
[TJ [2] it was shown that the 2 P°/ 2 and ^3/2 initial metastable states of the Se + ion 
both produce two dominant Rydberg resonance series. Our results for these resonances 
are illustrated in figure uj and M The 4s 2 4p 2 ( 1 D2) nd series converging to the X T>2 
series limit is indicated with solid triangles and the 4s 2 4p 2 ( 1 So) nd series converging 
to the ^o limit by the open triangles. Furthermore, two additional autoionizing 
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Figure 11. (Colour online) Effective quantum number u = n - /i and the 
linewidth V (meV) versus photon energy (eV) for the 4s 2 4p 2 ( 1 D2)nd resonances 
lying below the X D2 threshold of the residual Se 2 + ion. Experimental data 
(open red squares, ALS [HE]) and theoretical estimates (solid blue triangles, 
QB); (a) effective quantum numbers (u) for the 4s 2 4p 2 ( 1 D2)nd series originating 
from the initial Se + ( 2 D°,„) state; (b) effective quantum numbers {v) for 

the 4s 2 4p 2 ( 1 D2)nd series originating from the initial Se+( 2 D°,„) state; (c) 

4s 2 4p 2 ( 1 D2)nd, resonance linewidths (T) originating from the initial Se + ( 2 D° ,„) 

state and (d) 4s 2 4p 2 ( 1 D2)nd, resonance linewidths (r) originating from the initial 

)° ) 
J i/2> 



Se+( 2 D° ) state. 



3/2 and 2D 5 



metastable states, 



and 10 



Rydberg series were identified associated with the 2 Dg 

our theoretical results are shown by the solid triangles in figures 

The relationship between the principal quantum number n, the effective quantum 
number v and the quantum defect /x for an ion of effective charge Z is given by v = 
n — \x where the resonance position e r can be determined from Rydberg 's formula 

e r = £oo - Z 2 /v 2 , (2) 

and Eoo is the resonance series limit [H]. Tables 2-5 lists the principal quantum 
number n, resonance energy e r and quantum defect /i of the first few members obtained 
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Table 2. Principal quantum number (n), resonance energies E (eV) and quantum 
defect (/i) from experimental measurements of Se + [T] [2] compared with present 
theoretical estimates from the QB method. The Rydberg series originating from 
the 2 P2 /, metastable state of Se+ due to 4p — > nd transitions is tabulated. The 
uncertainties in the experimental energies are stated to be ±0.010 eV or less. 
The series limits are taken from the NIST 38 tabulations and from Esteves and 



co-workers lll[2l. 
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00 
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from our work for each of the prominent Rydberg series originating from the ^3/2 and 
2 P°/ 2 metastable states. The experimental results of Esteves and co-workers [TJ [2] are 
included in Tables 2- 5 for comparison purposes where it is seen our theoretical work 
is in satisfactory agreement with experiment for the resonance energies and quantum 
defects. Measurements taken at the Advanced Light Source synchrotron radiation 
facility indicated the presence of several prominent 4s 2 4p 2 ( 1 D2) nd resonances in this 
energy range originating from the 2 D£ , 2 and 2 Dg , 2 initial metastable states. Tables 6 
and 7 present results for the 4s 2 4p 2 ( 1 D2) nd dominant resonance series observed in the 
PI cross sections originating from the 2 Dg , 2 and 2 Dg , 2 metastable states respectively. 
Here again our theoretical work show satisfactory agreement with experiment for the 
resonance energies and quantum defects. 

Finally from figure [TT] theoretical estimates indicate many of the linewidths (r, 
determined from the QB method) for the 4s 2 4p 2 ( 1 D 2 ) nd series are less than a few meV, 
making their experimental determination challenging. Overall, given the complexity 
of this trans-iron element and the wealth of resonance structure in the corresponding 
PI cross sections, the agreement between theory and experiment is very satisfactory. 
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Table 3. Principal quantum number (n), resonance energies E (eV) and quantum 
defect (/i) from experimental measurements of Se + [T] [2] compared with present 
theoretical estimates from the QB method. The Rydberg series originating from 
the 2 P%/ 2 metastable state of Se+ due to 4p — > nd transitions is tabulated. The 
uncertainties of the experimental energies are stated to be ±0.010 eV. The series 
limits are taken from the NIST 38 tabulations and from Esteves and co-workers 
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4. Conclusions 



State-of-the-art theoretical methods were used to investigate the photoionization of 
Se + ions in the energy region of the ground-state ionization threshold (18 eV - 31 
eV). Given the complexity of this trans-iron element, overall, satisfactory agreement is 
found between theoretical and experimental results both on the photon-energy scale 
and on the absolute PI cross-section scale. The photoionization cross section exhibits 
a wealth of resonances. Prominent members of the Rydberg series are analyzed 
and compared with experiment. Resonance positions and quantum defects for the 
4s 2 4p 2 ( 1 D 2 )nd and 4s 2 4p 2 ( 1 So)nd series converging to their respective series limits are 
presented. The metastable state content of the Se + ion beam has been determined 
along with prominent autoionizing Rydberg resonances series seen in the experimental 
PI cross sections. A clear pattern of Rydberg resonance series is seen in the spectra 
of the photoionization cross sections disrupted by interloping resonances (see figures 
4 and 5). 

We point out that the strength of the present study is in the convergence between 
two different state-of-the-art theoretical methods and the comparison with available 
experimental data. Since no beam diagnostics were carried out in the experiments 
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Table 4. Principal quantum number (n), resonance energies E (eV) and quantum 
defect (/i) from experimental measurements of Se + [T] [2] compared with present 
theoretical estimates from the QB method. The Rydberg series originating from 
the 2 P° .. metastable state of Se+ due to 4p — > nd transitions is tabulated. The 
uncertainties of the experimental energies are stated to be ±0.010 eV. The series 
limits are taken from the NIST 38 tabulations and from Esteves and co-workers 
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[TJ[2] it is left to theory to infer the appropriate fractions. Such detailed comparisons, 
between theory and experiment, strengthens the validity of our results for use in 
astrophysical applications. 

The photoionization cross-sections from the present study are suitable to be 
included into state-of-the-art photoionization modelling codes Cloudy and XSTAR 
|19l I20| that are used to numerically simulate the thermal and ionization structure 
of ionized astrophysical nebulae. In combination with our determinations of the 
photoionization properties of other low-charge Se ions (to be addressed in subsequent 
publications), this will enable the abundance of Se to be accurately determined in 
astrophysical nebulae. Such information bears implications for the sites of origin and 
chemical evolution of Se and other trans-iron elements in the Universe, as well as 
for studies of the interior structure and nucleosynthesis of giant stars. These recent 
modifications to the Dirac-Atomic-R-matrix-Codes (DARC) have now made it feasible 
to systematically study other heavy complex systems of prime interest to astrophysics. 
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Table 5. Principal quantum number (n), resonance energies E (eV) and quantum 
defect (/i) from experimental measurements of Se + [T] [2] compared with present 
theoretical estimates from the QB method. The Rydberg series originating from 
the 2 P° .. metastable state of Se+ due to 4p — > nd transitions is tabulated. The 
uncertainties in the experimental energies are stated to be ±0.010 eV. The series 
limits are taken from the NIST 38 tabulations and from Esteves and co-workers 
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Table 6. Principal quantum number (n), resonance energies E (eV) and quantum 
defect (/i) from experimental measurements of Se + [T] [2] compared with present 
theoretical estimates from the QB method. The Rydberg series originating from 
the 2 D° ,„ metastable state of Se + due to 4p — > nd transitions is tabulated. The 
uncertainties in the experimental energies are stated to be ±0.010 eV. The series 
limits are taken from the NIST 138 tabulations and from Esteves and co-workers 
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Table 7. Principal quantum number (n), resonance energies E (eV) and quantum 
defect (/i) from experimental measurements of Se + [T] [2] compared with present 
theoretical estimates from the QB method. The Rydberg series originating from 
the 2 D°,„ metastable state of Se + due to 4p — > nd transitions is tabulated. The 
uncertainties in the experimental energies are stated to be ±0.010 eV. The series 
limits are taken from the NIST 38 tabulations and the work of Esteves and 
co-workers [Tll2]. 
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